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Abstract 

A thin water film evaporating from a cleaved mica substrate undergoes 
a first-order phase transition between two values of film thickness. During 
evaporation, the interface between the two phases develops a fingering 
instability similar to that observed in the Saffman- Taylor problem. The 
dynamics of the droplet interface is dictated by an infinite number of 
conserved quantities: all harmonic moments decay exponentially at the 
same rate. A typical scenario is the nucleation of a dry patch within the 
droplet domain. We construct solutions of this problem and analyze the 
toplogical transition occuring when the boundary of the dry patch meets 
the outer boundary. We show a duality between Laplacian growth and 
evaporation, and utilize it to explain the behaviour near the transition. 
We construct a family of problems for which evaporation and Laplacian 
growth are limiting cases and show that a necessary condition for a smooth 
topological transition, in this family, is that all boundaries share the same 
pressure. 
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1 Introduction 



The dynamics of evaporation and dewetting of thin liquid films arc important 
aspects of modern day technologies, such as in the fabrication of electronic chips, 
microfluidic devices and biosensors. These noncquilibrium processes depend on 
various details such as the temperature distribution in the sample, the nature of 
the interaction between the substrate and the fluid, and the substrate roughness 
[TJI2]- Depending on these details, a thin liquid film may, for instance, break 
into isolated droplets, evaporate uniformly, or exhibit a fingering instability 
similar to that observed in the Saffman- Taylor problem [3] where an inviscid 
fluid penetrates a viscous fluid. 

An experimental system showing this kind of behavior consists of a clean 
mica substrate covered by a thin water film [U |5J [5] . Under proper conditions 
of temperature and pressure, two films of different thickness (approximately 2 
and 12 nm) can coexist. This behavior is a consequence of the competing van 
der Waals and polar surface forces between the water and the substrate [7]. A 
first order phase transition between the film heights is induced by changing the 
vapor pressure. On the scale of nanometers, gravitational effects are negligible, 
and if the evaporation rate is sufficiently small one may neglect changes in 
the film thicknesses due to pressure gradients. This implies that the dynamics 
of evaporation, in this system is, effectively, that of the interface between the 
phases. 

An interesting phenomenon found in such a system is the nucleation of an 
'air' hole, such that an area of the thin film appears inside the thick film's area. 
This instability can then grow as the two boundaries of the thick phase move to 
meet one another. This can lead to breakup of the film into several components, 
or to the 'expulsion' of the air-hole from within the droplet. In topological terms, 
the 'droplet' (the thick phase) undergoes a topological transition from genus 
to genus 1, which can then be followed by other transitions: either back to genus 
0, or to a higher genus with some more complicated structure. 

Recently, the evolution of evaporating thin films was related to Hele-Shaw 
type [8| dynamics [9], adding it to a long list of physical systems whose dynam- 
ics can be approximated in similar fashion, for example the so-called oil-well 
problem |10j . the thin-sheet stamping problem and others. Such problems 
deal with the interface dynamics of two liquids with very different viscosities: 
A classic example is a droplet of a thin water film surrounded by air, which is 
a good description of the experimental setup described above. Many of these 
descriptions come under the general envelope of the idealized Laplacian growth 
problem. 

The connection to Hele-Shaw dynamics allows for the construction of an in- 
tegrable model for the interface dynamics of the evaporating films. Then, one 
can devise explicit analytical solutions with which to conveniently explore the 
properties of such systems. See for example |12| [5| 111! I13| among many others. 
In this paper we shall construct such explicit solutions describing the topolog- 
ical transitions described previously for the evaporating films. We will discuss 
the effect of the transition on the boundary conditions necessary for a complete 
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description of the problem. For example, adding an air-hole requires one to 
specify the pressure in the hole. This dramatically affects the behaviour of the 
interface as two boundaries approach one another and is closely related to previ- 
ous work regarding idealized Laplacian growth problems (see e.g. [14 ). Finally 
we shall characterize the transitions themselves and the dynamics close to such 
transitions. At this stage we will draw the connections with Laplacian growth 
problems. We will first demonstrate the duality of our system and Laplacian 
growth. Then we shall investigate how the stability of explicit solutions near 
the transitions depends on this duality. In the process we shall show how both 
Laplacian growth and thin film evaporation can be seen as special cases of a 
family of closely linked interface dynamics problems. This relationship becomes 
important when considering topological transitions. 

The outline of this paper is as follows: In section [2] we introduce the the- 
oretical framework necessary to describe the phenomenon we are discussing. 
In section [3] we describe how to construct instructive solutions using methods 
borrowed from the idealized Laplacian growth problem. In section [4] we shall 
bring two concrete examples of such solutions. In section|5]we shall examine the 
topological transitions in detail. Finally, in section [6] we will draw connections 
between transitions in evaporation problems and in Laplacian growth problems. 
We shall use these connections to explain some of the features discussed in 
section [5] 



2 Theoretical model 

The first step to modeling our system is to divide it into two regions, according 
to the two phases described above. Since both phases have very small thickness, 
the dynamics are controlled by viscosity. By making a no-slip requirement along 
the substrate-liquid interface, and a no-shear requirement along the liquid-air 
interface, we can average over the vertical direction, arriving at the D'arcy law 
for the liquid velocity v: 

h 2 

v = T VP, (1) 

where h is the film height, v is the viscosity, and P = P(r) is the average 
pressure at a point r = (x,y). The liquid dynamics in the thin phase may 
be ignored, since as ([TJ shows the velocity there is approximately zero. We 
therefore call the thicker phase a 'droplet' or water region and the thin phase 
a 'dry' or air region. Henceforth we denote the droplet region by the symbol 
D. Next we assume that the film height in the droplet region is approximately 
fixed, allowing us to neglect variations in the two-dimensional liquid density p. 
Uniform evaporation then implies 

pV ■ v = -k, (2) 

where k is the rate of evaporation. Combining and Q we see that the 
pressure obeys Poisson's equation in the droplet. For the time being, let us 
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ignore the possible nucleation of an air-hole in the droplet. Then without loss 
of generality we may take the pressure in the dry domain to be constant and 
equal to zero. Thus: 

V 2 P = f3 inside D, 
P = outside, ^ ' 

where j3 = 2>kv / ph 2 . From now on we shall assume units where it is dimension- 
less, along with both distance and time. A last simplification is to assume that 
surface tension can be neglected along the interface between regions so that the 
boundary condition on ^ is: 

P(r <= dD) = 0. (4) 

Here (and henceforth) dD denotes the droplet boundary (for both simply- and 
doubly-connected droplets) . The last assumption Q has the effect of endowing 
our system with an infinite number of first integrals of motion [T4l ITT] . Indeed, 
let h be any harmonic function on D. Then we may define the following quantity: 

t h = [ h(r)d 2 r. (5) 

JD 

We find the time dependence of th' 



d t h 
dt 



d -fhdl 
,, D On 

— h dl + 

dD on Jg D 



on 



[ (-hV 2 P + P\J 2 h) d 2 
Jd 

-j3h d 2 r = -Pt h 



(G) 



In this derivation we have used Green's theorem along with (Tl and Q. d/dn 
is the derivative perpendicular to the boundary. To utilize ( 6 1 , we define the 
so-called harmonic moments: 



t k = z dx dy, 
7T Jd 

t k (t)=t k (0)exp(-pt). 



(7) 



The time dependence of these moments can be used, in principal, to completely 
determine the interface dynamics. Let us clarify this through an example, by 
solving the problem of a circular droplet of radius L with initial radius Lq. 
Solving ([3|+(|4| in polar coordinates (r, 6) we obtain P(r,t) — /3/A[r 2 — L(t) 2 }. 
In order to find L(t) we note that tk = for any k > while to — L(t) 2 = 
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Lpexpf— fit). Thus L(t) = L exp(— fit/2) solves the problem. Now let us see 
what happens once we allow the creation of an air-hole in the droplet, which 
is the case this paper treats extensively. Let us use our previous example to 
illustrate. Assume that an air-hole of radius Iq nucleates in the center of our 
circular droplet at t = 0. The solution of ([3]), Q now changes to 

P(r) = A+ ^r 2 + -logr (8) 

where A = ~fil 2 /A -B\ogl/A,B = B (l 2 - L 2 ) /\og(L/l). We see that adding 
the air-hole, i.e. changing the the droplet from being simply-connected to being 
doubly- connected has added a degree of freedom to the system, no longer de- 
termined by the harmonic moments Q alone. This is clearly a general aspect 
of such a change, not just an artifact of our example. We must constrain this 
freedom by adding an additional moment, which must express this topological 
change. This can be, for example, the logarithmic moment [14] 

U = - I log |r| dxdy. (9) 
71" Jd 

ti decays exponentially like the other moments, since ^o<?|r| is harmonic as 
long as ^ D (as we shall assume henceforth for simplicity). In the case 
of our evaporating droplet, now in the shape of an annulus, this furnishes 
ti = L 2 log L — I 2 log I — L 7 l . Adding this to t — L 2 — I 2 gives two non- 
linear equations for I, L, solving the problem. As expected, this solution depicts 
an annulus whose inner and outer boundaries approach each other. The solution 
exists for any finite t and approaches a circle. One can find the circle's final 
radius by taking the limit I —> L, which yields Lit — > oo) = exp[f;(0)/£o(0)]. Of 
course, for such a simple geometry one can solve the problem by direct applica- 
tion of |8]) to DArcy's law Q, without recourse to moment equations |15j . 

It should be noted, that when we claimed that adding an air-hole doesn't 
change the time dependence of the moments ([6]), we implicitly used our as- 
sumption that the pressure is equal on both the inner and outer interfaces of 
the droplet, that is 5P = between boundaries. Although this certainly applies 
to the physical system we have in mind, one can readily think of cases where 
this is not so. For example, we can envision a system where (|3| applies yet the 
pressure is zero only on one boundary. This could be a case where a droplet is 
compressed between two solid sheets at an exponential rate. An alternative is 
to envision that we take the limit fi — > but allow the pressure in the hole to 
be different from outside. This would be relevant for a Hele-Shaw cell, where 
we pierce the top sheet with a small hole and use that to enforce some pressure 
difference between boundaries. In this case, re-using the arguments in dfjl, the 
time dependence of the moments becomes: 

^ = / (- hV 2 P + PV 2 h) d 2 r- 5 ^i ^dl 
at Jd tt J dDout an 

sp r 

= -fit h / V 2 /i d 2 r. (10) 

77 JD out 
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Here dD out and D out are respectively the outer boundary of the droplet and the 
area it encloses. We have assumed without loss of generality that the pressure 
on the inner boundary is zero, while it is SP on the outer boundary. Since z 
are analytic in D out , the t^'s time dependence remains (|6]). but now 

f— »-MP. (11) 

since V 2 log(r) = 27r<5 (2) (r). 

Let us briefly summarize this section. We have simplified the problem of an 
evaporating droplet to solving Poisson's equation for the pressure, with P = on 
the droplet boundary, for a moving boundary. This formulation leads to simple 
laws for the harmonic moments of the droplet, and these can be used to extract 
the time dependence of the boundary. A topological transition in the droplet, 
for example the appearance of an air-hole, adds a degree of freedom. This 
freedom has the physical interpretation of the ability to determine the pressure 
in the air-hole by various external means. It requires adding a constraint in the 
form of the logarithmic moment, whose time dependence encapsulates our new 
physical information on the system. Once this moment is known, we can again 
solve the given problems. 

For a more extended discussion of these points, and those in the following 
section, we refer the reader to the large body of work on Laplacian growth, for 
example |16llT7] or the review [IB], as well as the works cited above. These also 
supply the connection between our problem and many other fields, which in the 
interest of brevity we do not address here. 



3 Method for constructing explicit solutions 

Unfortunately, finding the interface dynamics from the infinite number of mo- 
ments tk (and ti) is in general a formidable problem. It can be simplified greatly 
by assuming that the area of the droplet D can be obtained as a conformal map- 
ping from some simple region, with known singularities and branch cuts. For 
example we can choose a mapping so that our droplet domain D (for now with- 
out a hole) is obtained as a mapping z(w) from the unit circle w < 1. In this 
way the moment constraints translate into equations for the parameters of z(w). 
To see this we introduce the Schwarz function of D, defined as a function that is 
analytic in a strip about the boundary, and on the boundary obeys the relation 
z = S(z). For a mapping z(w) as defined above, we have 



S(z) 



1 



w(z) 



(12) 



Here z(w) = z(w) is the conjugate to z(w), and the 1/w term appears since on 
boundary of the unit circle w = 1/w, and we continue analytically from there. 
w(z) is the inverse map of z(w), which in general is multivalued, so that S(z) is 
analytic on a (multi-sheeted) Ricmann sphere . If z(w) has poles of order rij at 
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Figure 1: (color online) Mapping from the unit circle to a droplet. The left 
panel shows a mapping z(w) that is conformal in the unit circle, having all its 
poles and branch points in the domain \z\ > 1. The contours represent lines of 
= Const while the color shading denotes the phase Arg z{w). The right 
panel depicts how the unit circumference is mapped to the droplet' domain. 
Embedded in the droplet are the poles of the Schwarz function. Note: the left 
panel shows the mapping for the largest droplet contour. 



points wj, by definition in the region w > 1, then S(z) must have singularities 
at points qj G D inside the droplet, and obeying: 



9j 



Wj 



to 



lim Sz nj z 



1 



w (qj)Sz 



(13) 



This allows us to rewrite ^ using Green's theorem as: 



z k S{z)dz 



27" J 3D 

J2to 



k-\-l — 7ij 



(14) 



It is clear that in order to obey ([6]), we must require that the poles qj be constant 
in time, while the weights fj,j all decay exponentially in time at the rate j3: 



qj = const, 



(15) 



Figure [T] depicts an example of such a mapping. 

In order to incorporate an air-hole, we must find a mapping from a different 
domain. The most obvious would be to require that D be obtained as a mapping 
z(u) from some annulus I < \u\ < L, such that each boundary of D maps to 
a boundary of the annulus [19 . This adds a degree of freedom, namely the 
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Figure 2: (color online) Mapping from a torus to a droplet with an air-hole. 
The left panel shows a mapping z(u) that is analytic in the right half of the 
rectangle, having all its poles and branch points in the domain u < 0. The 
contours represent lines of = Const while the color shading denotes the 

phase Arg z(u). The mapping is elliptic, i.e. one should consider the top/bottom 
and left/right boundaries to be connected and forming a torus. The verticle 
lines 5Ru = 0, 5ft u = a x are mapped to the droplet boundaries. Embedded in the 
droplet are the poles of the Schwarz function. 



modulus of the annulus [20J, which is then constrained by £/. In practice, we 
found it easier to work instead with a rectangle, so that z(u) is an elliptic 
mapping from a rectangle of half- widths a x ,a y . This means that S(z) = z(—u), 
so that S(z) is analytic on a torus of modulus r = a y /a x . See figure [2] for an 
example of such a mapping. Now, to rewrite tf. via Green's theorem we must 
introduce a cut in D, from some point on the inner boundary zq to some point 
on the outer boundary Z\ (see figure |3|. It is easy to see that ( [l4| is not changed 
by this addition. However, applying Green's theorem to t\ requires care, as D 
now winds around a singularity of log(z). We write t\ = 5ft ^ J D logz d 2 z, and 
then using Green's theorem we have: 

2 1 I 

log z d z = </> log zzdz 

2m J 

1 1 log zS(z)dz + i [ f +( ) [z-S(z)]hgzdz 

(16) 



2ni Jq D 2iri 



where the final term describes integrals along both branches of the cut (see 
figure |3j. The first term on the right clearly depends only on fij, qj. To evaluate 
the second term we note that since log z = log \z\ +i arg z, the two integrals with 
the integrand (. . .) log \z\ exactly cancel, and with proper choice of the cut, one 
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Figure 3: Integration path for the moments 



integral is on the branch arg z — and the other is on the branch arg z = 2it. 
Thus: 



j 



S{z)]dz 



(17) 



Here Rj are simply the appropriate residues of log 2. The first term decays 
exponentially. The second term must decay at the same rate, and this provides 
the constraint on the torus modulus r that is necessary to complete our system 
of equations. In fact, we may define the so-called 'modified Schwarz potential' 

ED 



Q(z) 



S{z')dz' - / S{z')dz' 



(18) 



named thus because it gives the contours z — S(z) as its saddle manifolds. 
Then, if we define a potential mismatcl^ 



SQ = 5R{Q(zi) - Q(z a )} = f 

J Z, 



S(z)]dz 



(M 



*0 



S(z)dz 



we may use the following constraint instead of (111: 

-P8Q- 



dSQ 

~dT 



SP 



(19) 



(20) 



since the second term on the right-hand side of (17 1 is precisely —2SQ. Fur- 
thermore, let us recall that S(z) is analytic on a Riemann surface of genus 1. It 
can be shown [23] that for a closed cycle T on the surface, going from zq to Z\ 



lr The use of the name 'modified Schwarz potential' is perhaps slightly misleading. Usually 
this function is chosen so that both it and its derivative are zero on the boundaries (i.e. 
SQ = 0). However, the quantity we use captures the essential relation of the modified Schwarz 
potential to the Schwarz function, as defined in 1211 . so we use it here. 
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in the droplet and then via an arbitrary route through the air region back from 
z\ to zq, one has: 



SQ 



S{z)dz. 

z{—u)z' \u)du, 

(r) 



(21) 



and we know that u(T) begins on the line 3?it — ~a x and goes to = +a x . 
Since S(z) is analytic on the surface, we may choose any loop homotopic to the 
one we chose, and we see that our choice of Zq,Z\ and Y was arbitrary. Note 
that the term ■potential mismatch is apt in the sense that two boundaries cannot 
meet unless SQ = 0, since SQ is defined in ( 19 1 as an integral from one boundary 
to the other. Furthermore, if the pressure difference is zero then the mismatch 
decays as the boundaries approach. However, from (20) we can immediately 



conclude that if there is a potential mismatch at the droplet 's initial state, the 
two boundaries will not meet in any finite time. 



4 Examples of explicit solutions 

Using the method we described, one can construct solutions with features of 
interest. We shall present two examples of such solutions. The first can be 
seen as describing the nucleation of a small air-hole in the droplet, which then 
expands to meet the outer boundary, finally combining so that the final droplet is 
again simply connected. To this end we seek a Schwarz function whose behaviour 
in the droplet is of the form 

S(z)=t^T + --- ( 22 ) 
qj G D, fii> 0. 

where (...) denotes some analytic function. This form has been used previously 
in a number of interesting cases of doubly connected regions in various problems, 
e.g. |24l [231 [231 127] , The need for three poles can be explained qualitatively 
as follows: If S(z) has one pole only, then the droplet it describes is highly 
symmetric, and there is no preferred spot for hole nucleation. The second pole 
therefore breaks this symmetry, and the third provides for the actual appearance 
of an air-hole. It is also possible to show formally |281 127] that three poles is 
the minimum necessary for air-hole nucleation. To obtain such a mapping, we 
define an elliptic mapping with half periods a x ,a y : 

2 

z(u) = (XjC(u — Uj) + 7; —a x < tftuj < 0, — a y < ^suj < a y (23) 

3=0 

where ((z) is the Weierstrass Zeta function, which is quasi-periodic. With prop- 
erly chosen parameters, the mapping is conformal for u in the right half of the 
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Figure 4: Nucleation of air-hole in a water droplet. The left panel shows the 
outline of the droplet defined by the mapping ( 23 ) . Although at that magnifi- 
cation it can't be seen, there is a tiny air-hole on the y = axis near x = 1.01. 
This can be seen in the magnified picture in the right panel. The contours in 
this figure denote equal-time steps in the evolution of the droplet domain. There 
are such contours in the left figure as well, but they are so closely spaced as to 
be indistinguishable. This shows how the short-scale dynamics of the air-hole 
evolution are much faster then that of the droplet as a whole. 



fundamental region of the mapping, thus defining the droplet domain. To ensure 
periodicity we demand 



o., 



(24) 



and along with (13), (17 1, ( |20| and ( |21) this provides exactly 8 constraints 
for the seven parameters in the mapping plus the modulus r = a v /a x . These 
nonlinear equations can be solved numerically. An example of such a solution 
appears in Figure if 



Our second example can be seen as describing how the nucleation of an 
air-hole ends up in breaking the droplet into separate fragments. We choose a 
mapping of the form [121 |2U] : 



z(u) = ^<x, [logcr (it - Uji) - log cr(u - u j2 )] + 7; 
3=0 

— a x < $iujk < 0, — a y < ^sujk < a v 

2 

S(z) = Y^ [log(z - qji) - log(z - q j2 )} + ... 

3=0 



(25) 



(26) 



2 The parameters used here as initial conditions were, up to 5 digits of accu- 
racy: (u , Ul ,u 2 ) = (-4.60281,-4.10607,-3.21041), (a , ai ) = (0.99312, -1.28713), 7 = 



0.99975, a x 



5.0, t = 1.2. 
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Figure 5: Breakup of a droplet as a result of air-hole nucleation. The figure 
denotes the evaporation of a droplet defined by a mapping of the form (251. 
The contours denote equal-time steps in the evolution of the droplet domain. 



where a(z) is the Weierstrass Sigma function. It can be verified that the re- 
quirements ( |15[ ) remain unchanged for a mapping of this form. The periodicity 
of z{u) and the quasi-periodic behaviour of the Weierstrass Sigma function lead 
to the requirement, 



Y^otjiuji - u j2 ) = 0, 



(27) 



which along with the nine equations (151 plus (20 1 fully determine the eleven 
parameters in ( 25 ) . Figure [5] denotes the evolution of such a droplet [^] For the 
examples shown in this figure, we have chosen a highly symmetric configuration, 
since otherwise it becomes very hard to numerically solve the large number of 
resultant nonlinear equations. 



5 Characterization of transitions 

In the previous section, we showed how explicit solutions exhibit features appear- 
ing in evaporating droplets. Both these examples showed how an evaporating 
droplet eventually undergoes a topological transition from genus 1 to genus 0. 
This contrasts sharply with the example of an evaporating annulus from section 
[2j where the droplet remains in the form of an annulus at any finite time - a 

3 The parameters used to start this solution were, to 6 digits of accuracy, Ujk = —4.4 + 
S k2 0.193 + i 14{3 3 ~ 1) ,aj = exp - 1),7 = 1.0, a x = 5.0, r = 1.2. 
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AQ<0 / \ AQ=0 / AQ>0 




Figure 6: (color online) Cusp formation due to pertubation of initial conditions. 
The middle panel is a depiction of the final frame in figure [4] The right and left 
panels depict the evolution of the droplet when SQ ~ ±10~ 7 . For comparison, 
all other parameters in this configuration take on values between an order of 
1CP 3 to an order of 1. 



clearly unphysical solution resulting from symmetry. One mathematical reason 
for this can be found immediately. For both examples in section [4] the initial 
conditions were chosen so that the potential mismatch SQ = initially (and 
thus forever). For the annulus, we can find the mismatch by noting that the 
Schwarz function for an evaporating annulus has, near the boundaries, the form 
R 2 /z, where R = l,L on the inner and outer boundaries respectively. Using 
this, one can easily identify that for the annulus ti — —SQ. This value is clearly 
non-zero as can be seen from the calculation in section [2] However, this dif- 
ference is not the whole story. Figure [6] depicts the result of pertubating the 
initial conditions of the mapping shown in figure [3] so that SQ ^ but is still 
very small compared to all other parameters in the problem. As can be seen 
in the figure, the result is a formation of a cusp on one boundary (leading to a 
breakdown of the solution) with the other boundary remaining regular. We note 
that a similar behaviour to the one shown in figure [6] was found in a different 
context in |14J . In addition, for this specific example, the velocity at the regular 
boundary goes to zero as the cusp forms. This behaviour is very similar to what 
one would expect were there actually a pressure difference on the boundaries. 
A simple explanation for this behaviour can be found by noting that as can 



be seen from (18 1, the function f3Q(z) itself solves Poisson's equation in the 
droplet domain. It cannot however be identified with the pressure, since it has 
singular points associated with the singularities of S(z). Nevertheless, near the 
transition area, far from the singular points, Q itself behaves very nearly like a 
pressure field itself, with the important addition that VQ = on the boundary, 
i.e. the boundaries are saddle manifolds of Q. When SQ = 0, then there is no 
problem for these saddle manifolds to remain stable until the actual transition 
when the boundaries meet. However, if SQ ^ 0, then one of two options must 
occur: either the dynamics slow down to ensure that the saddle points remain 
stable, or if the interface velocity is too high then at some point one of the saddle 
manifolds, associated with the higher value of Q, must become unstable, leading 
to a breakdown of the solution. Thus Q acts as a quasi-pressure field itself, and 
the difference in this quasi-pressure causes the breakdown of the solution. 
We shall now investigate this pressure analogy further. 
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6 Duality of Laplacian growth and evaporation 



In this section we shall give the modified Schwarz potential Q an interpretation 
as a real pressure field. We shall do this by utilizing a duality between thin film 
evaporation and the idealized Hele-Shaw cell, or Laplacian growth, problem. 
By duality we mean that for a given evaporation problem with some known 
solution, we can construct an analogous Laplacian growth problem leading to 
a different pressure field but with the same solution for the interface dynam- 
ics. This relation has also been studied in the context of the 'inverse balayage' 
problem in potential theory and quadrature domains (see e.g. [30] and |31|). 
The pressure field in this dual problem will naturally contain Q as one of its 
elements. We will then generalize this duality to a family of analogous prob- 
lems, all featuring the exact same interface dynamics. This generalization will 
demonstrate how a potential mismatch, i.e. 5Q =^ 0, in an evaporation problem, 
translates to a pressure difference in the analogous problems, without affecting 
the interface dynamics, and explains the behaviour shown in figure [6j 

Let us briefly cover the definition of the idealized Hele-Shaw cell evolution 
problem. Suppose that instead of evaporation, the droplet dynamics are gov- 
erned by the existence of pumps at discrete points qj , such that in the vicinity of 
these points the pressure diverges as P ~ log \ z — In this case, the pressure 
obeys the following relation: 

V 2 P = Y,Rj{t)S (2) {z-<l3) (28) 

3 

where Rj(t) denote the rate of pumping at the points qj. Positive Rj are 
equivalent to liquid extraction, and negative Rj to injection. In the case of 
a simply connected droplet, it is possible to reduce our evaporation problem 
to Laplacian growth, simply by placing sinks with rates Rj{t) = 0Trfij(t) at 
the poles of the Schwarz function. This can be verified by directly calculating 
the harmonic moments by methods analogous to those used in ([6]) and (14 1 



However, in the case of a multi-connected droplet, in order to obtain the same 
evolution for the potential mismatch, and thus the same dynamics, one must 
apply a pressure difference of —f36Q(t) on the boundary. This can be seen 



from (20 1, by taking the limit (3 — > 0. Thus we see that the dynamics of an 
evaporating droplet with potential mismatch SQ are the same as a Hele-Shaw 
droplet with a pressure mismatch —(35Q. In fact, by using the correspondence 



of the singularities of S and Q, as is evident from (18), we can readily identify 



that the solution to ( 28 ) with the relevant boundary is simply 

P = P-PQ. (29) 

With this, the fact that the dynamics are the same for this droplet is trivially 
obtained from Q's saddle-manifold behaviour: the interface velocity is precisely 
v oc —VP = —VP. For SQ ^ the pressure difference is non-zero and so 
we know that the boundaries cannot meet. In this case it is less clear how 
the pressure difference causes the cusp formation, except that one expects that 



14 



the area of higher pressure will actually be more regular (as opposed to the 
evaporation problem) as is indeed the case. We still don't have a clear picture 
of why the regular boundary should slow down. 

We can take this point further by trying to construct more problems featur- 
ing the same dynamics. To do so, we can envision a situation where we have 
both evaporation and pumping as evolution mechanisms. Let us concentrate on 
the case where we are moving liquid via pumps, at an exponentially decreas- 
ing rate, while evaporating at a different rate, so that the pressure obeys the 
condition 



V 2 P = (3 + a 



S^(z-q j} 



(30) 



so that the droplet is evaporating at a rate j3 + a, but the a evaporation is 
being exactly canceled by injection (or extraction). Clearly a — reduces to 
our original problem, and a — —ft reduces to Laplacian growth. Now, for a 
harmonic function h we have, by adapting (10) to our new system 

SP 



d t h 
dt 



~hV 2 P + PV 2 h) d 2 r- 



-h{(3 + a 



P\7 2 h d 2 r 



dh 
On 



dl 



S^(z 



1i) 



d 2 r- 



(31) 



SP 

7T 



V 2 h d 2 r 



For the moments tk, i.e. h — z k /-k, the last two terms on the RHS fall and the 
first simply evaluates to — fit^, from (14 1. However, for ti, that is h — tt^ 1 ^. log z, 
we have: 

^ = - (/3 + a )ti + a M j (*) log 9, - 25 P. (32) 



Using (17) we can rewrite this as: 

-(J3 + a)ti + a(ti-2SQ)-2SP, 



dti 
dt 



(33) 



and so in order to get the same dynamics as the evaporating droplet we need: 

aSQ(t) = 5P(t) (34) 

So, for any a / we see that applying a very specific pressure difference will 
yield the same dynamics. Note, that in this case the pressure has to be deter- 
mined from Q and not vice- versa. (This is in close analogy to the case discussed 
in |23|.) In addition, we see that any a > is equivalent to having a ■positive 
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pressure difference which is exactly the behaviour that explains the shapes in 
figure [6] For any a < we have the equivalent of negative pressure difference. 
In that sense, the question of whether the pumps are extracting or injecting 
liquids appears to determine the mechanism of breakdown. Finally, it is again 



clear that the solution to (30 1 is precisely P = P + aQ. In this sense Q repre- 
sents the pressure field resulting from evaporation along with a compensating 
injection. Our original analogy indeed manifests itself in a concrete manner. 

We should emphasize a couple of points. First of all, we have not shown that 
a pressure mismatch 5Q ^ necessarily leads to one of the boundaries coming 
to a halt. The solution to a given problem could well break down long before 
one boundary comes to a halt. In other words, a breakdown will occur when a 
saddle contour of Q becomes unstable, but that doesn't mean that the interface 
velocity will have reached zero by then. We have also not shown that SQ = 
must lead to a meeting of boundaries even if the only important dynamics are 
occuring near the transition point. Perhaps cusps appear before the meeting. 
We have only shown that the mismatch leads to a bound on the time for which 
a solution exists, and provided a mechanism for the breakdown. The actual 
exact relation between a solution's breakdown and the evolution of Q remains 
a matter for further study. 
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